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Brief, Popular Summary of the Paper: 

Atmospheric computer models have been used for about the last 30 years to help 
researchers understand the response of stratospheric ozone to various natural and human-made 
variations. Two-dimensional (latitude versus altitude) models have been developed in this time 
period and include a representation of transport and most of the chemical processes important in 
driving global ozone changes. Both the atmospheric transport and chemistry are important in the 
computation of ozone levels, thus calculated ozone variations from human- induced changes are 
dependent on both the transport and chemistry representation in the individual model. 

The three two-dimensional (2-D) models AER, LLNL, and GSFC were developed at 
Atmospheric and Environmental Research in Cambridge, MA, Lawrence Livermore National 
Laboratory in Livermore, CA, and Goddard Space Flight Center in Greenbelt, MD, respectively. 
These models have been used extensively since the mid-1980’s to assess the influence of various 
human-made impacts on stratospheric ozone. 

Until this study it has been difficult to understand the differences in results from these 
three models since each model was developed separately by different groups of investigators and 
were dependent on both the 'ransport and chemistry representation. The AER team was able to 
develop a common framework a few years ago that included the most relevant parts of the LLNL 
and GSFC models. With this common framework the three models could be compared much 
more quantitatively and the i ransport and chemistry portions of each model could be separated 

and compared individually with one another. 

Differences in the model simulations of long-lived gases pointed out the transport 
differences among the models and the variations in the simulations of families of more reactive 
gases pointed out the chemistry differences. Simulations of the influence of the proposed High 
Speed Civil Transport (HSCT) planes indicated huge transport and cold polar process driven 
differences among the models. 
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Abstract. 

A model intercomparison among the Atmospheric and Environmental Re- 
search (AER) 2-D model, the Goddard Space Flight Center (GSFC) 2-D model, 
and the Lawrence Livermore National Laboratory (LLNL) 2-D model allows 
us to separate differences due to model transport from those due to the model’s 
chemical formulation. This is accomplished by constructing two hybrid mod- 
els incorporating the transport parameters of the GSFC and LLNL models 
within the AER model framework. By comparing the results from the na- 
tive models (AER and e.g. GSFC) with those from the hybrid model (e.g. 

AER chemistry with GSFC t ansport), differences due to chemistry and trans- 
port can be identified. For tb 1 analysis, we examined an inert tracer whose 
emission pattern is based on emission from a High Speed Civil Transport (HSCT) 
fleet; distributions of trace species in the 2015 atmosphere; and the response 
of stratospheric ozone to an HSCT fleet,. Differences in NO y in the upper strato- 
sphere are found between models with identical transport, implying differ- 
ent model representations of atmospheric chemical processes. The response 
of 0 3 concentration to HSCT aircraft emissions differs in the models from 
both transport-dominated differences in the HSCT-induced perturbations 
of H 2 0 and NO y as well as from differences in the model representations of 
0 3 chemical processes. The model formulations of cold polar processes are 
found to be the most significant factor in creating large differences in the cal- 


culated ozone perturbations. 
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1. Introduction 

Studies have shown that transport is a major uncertainty in 2-D and 3-D modeling of the 
atmosphere [ Jackman et al., 1991; Kinnison et al., 1994: Douglass et ai, 1999]. Models 
also differ in their treatment of chemistry. Though most, models employ the chemical 
reaction rates from JPL [DeMore et al., 1997; Sander et al., 2000], differences exist in 
family groupings, diurnal avei aging technique, and treatment, of heterogeneous chemistry. 
Calculations of ozone perturl -ations by different models show differences which are not 
easy to interpret. How much of these differences are due to differences in transport? How 
much to differences in chemistry? Better understanding of intermodel differences would 
help in placing uncertainty estimates around model predictions of future ozone changes. 

Model intercomparisons have typically involved comparisons of spatial and temporal 
distributions of trace gases. Information on differences in chemistry is obtained by com- 
paring radical concentrations calculated using specified long-lived species concentrations 
(see e.g. the chemistry intercomparison exercises in M&M I [Prather and Remsberg, 1993] 
and M&M II [ Park et al., 199!*]). Intercomparisons of model dynamics use idealized tracer 
species calculated with prescribed production and loss rates. Tracer concentrations are 
then compared, or quantities such as age of air are derived from the tracer mixing ratios 
[Hall et ai, 1999; Park et al, 1999]. 

This paper will employ a different approach. This study is based on three two- 
dimensional (2-D) models: the Atmospheric and Environmental Research (AER) 2-D 
model, the Goddard Space Flight Center (GSFC) 2-D model, and the Lawrence Liv- 
ermore National Laboratory LLNL) 2-D model. The transport parameters, and to the 
extent possible, the numerical scheme and grid resolutions of the LLNL and GSFC models 
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have been implemented within the AER model. Using this approach, chemical differences 

for both long- and short-lived species can be evaluated by comparing two models with 
similar transport but differem chemical formulations. The effect of transport on model- 
calculated perturbations can be evaluated by comparing results from the AER model 
with different transport formulations but identical chemistry. This modular approach to 
model intercomparison is beii g used in the 3-D Global Modeling Initiative (GMI) model 
[Douglass et al, 1999; Considme et al, 2000; Rotman et al, 2001; Kinnison et al, 2001]. 

The following section will describe the dynamical and chemical formulations of the 
three models. Section three will describe the intercomparison approach used. Section 
four will present an intercom parison of a background atmosphere calculation for 2015 
conditions. Section five will present an intercoinparison of calculated perturbations to 
H20, N(\, and ozone due to a future fleet of stratospheric aircraft. Section six includes 
discussion and conclusions. It should be noted that the results presented in the paper were 
generated by the models around 1997 in preparation for the M&M II [Park et al., 1999] 
and IPCC [1999] reports. Each model has continued its modifications and refinements 
since then. Although many < etails of the results presented in this paper may well have 
been superseded by later wo)k, this in-depth comparison of three 2-D models provides 
insights into many intermodel differences. This paper also presents a demonstration of 
a new capability for comparing different 2-D model results, in which one of the models 
provides the framework for inclusion of the advection and diffusion of other models. 

2. Model Descriptions 

The three models considered here differ substantially in model formulation and details 
of chemistry and dynamics. All three models use log pressure as the vertical coordinate, 
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but vertical resolution varies f om 1.2 km in the AER model to 1.5 km in the LLNL model 
to 2.0 km in the GSFC mod* A The GSFC model uses 10 degree horizontal resolution, 
the AER model 9.5 degree hoi izontal resolution, and the LLNL model 5 degree horizontal 
resolution. The top of the GS '0 model extends to 90 km, the LLNL model to 80 km, and 
the AER model only to 60 km All three models use observed climatological temperatures 
taken from NCEP analyses [ I'alnay et al., 1996]. 

2.1. Transport Parametei s 

The transport coefficients o the LLNL model [Kinnison et al., 1994; Li et al., 1995] are 
computed using model-calculated ozone and observed temperatures. Mechanical forcing 
from planetary waves breaking (represented explicitly by two waves), gravity wave break- 
ing, and Rayleigh friction are also employed. While this calculation is done interactively 
in the original LLNL model, the transport parameters from a previous calculation have 
been saved and used repeatedly in the LLNL results arid the hybrid model calculations 
presented here. The horizon' al diffusion K yy is calculated as the ratio of the E-P flux 
divergence to the gradient of potential vorticity [ Garcia . , 1991], with a minimum value set 
at lxlO 8 cm 2 /s in the stratosphere. Tropospheric K yy is set to lxlO 10 cm 2 /s. A gravity 
wave breaking parameterizati >n is used to calculate K zz in the stratosphere, with values 
ranging from 2xl0 3 cm 2 /s to over lxlO 5 cm 2 /s. K zz in the troposphere is set to 4xl0 4 
cm 2 /s. Figure 1 shows K yy and vertical velocity from the LLNL model for June and De- 
cember. Note the large gradients in the K yy fields, especially in December where values 
of the equatorial K yy exceed 40xl() 10 cm 2 /s near 40 km. 

The GSFC model (see J adman et al. [1996]; Fleming et al. [1999]) calculates trans- 
port, parameters from observed climatological values of temperature, H 2 O, zonal wind 
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and ozone. Mechanical forcin ' from six planetary waves (constructed from the observed 

temperature field) and effect:, of gravity wave breaking provide the wave driving. The 
diabatic heating rates are computed following Rosenfield et al. [1994], and latent heating 
in the troposphere is taken from Newell et al. [1974], The streamfunction is obtained by 
solving an elliptic equation obtained by combining the zonal-mean momentum arid energy 
equations (see e.g., Garcia an 1 Solomon [1983]). K yy is computed using a similar theoret- 
ical basis as that used by the LLNL model but employing observed meteorological data 
and following the approach ol Randel and Garcia [1994]. K yy values range from lxlO 8 to 
5xl0 10 cm 2 /s, with values computed in both troposphere and stratosphere. K zz values in 
the troposphere and lower sti atosphere are based on the vertical temperature gradients. 
In the upper stratosphere and mesosphere, K zz values are obtained from a gravity wave 
parameterization. Stratospheric values of K zz range from lxlO 2 cm 2 /s in the lowermost 
stratosphere to 5xl0 5 cm 2 /s in the upper stratosphere. Tropospheric values of K, z range 
from lxlO 4 cm 2 /s at the tropopause to 5xl0 5 cm 2 /s at the surface in the tropics and 
2.5x10 s cm 2 /s at the surface tear the poles. Figure 2 shows K yy and vertical velocity for 
June and December from the GSFC model. 

The AER model uses transport parameters obtained in a more ad hoc manner than the 
GSFC or LLNL models. He* ting rates, based loosely on Dopplick [1979], are scaled by 
temperature lapse rate and used to obtain the vertical velocity. The vertical velocity field 
is integrated to obtain a stre unfunction, with adjustments made to assure mass conser- 
vation (see Ko et al. [1985]). Diffusion rates are specified independently of the advective 
circulation, with K yy values in the lower stratosphere based on studies of exchange time 
scales between the tropics and midlatitudes [Shia et al., 1998]. Values chosen are from 0.7- 
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1.3xl0 9 cm 2 /s in the tropics end from 3-10xl0 9 crn 2 /s elsewhere. K zz values are constant 

at lxlO 5 cm 2 /s in the troposphere, lxlO 3 crn 2 /s in the lower and middle stratosphere, 
and lxlO 4 cm 2 /s in the upp< r stratosphere. K yz values are obtained by projecting the 
K yy values from isentropic surfaces to pressure surfaces. There has been no attempt to 
maintain consistency between the advective and diffusive components of transport within 
the AER model. Figure 3 shows K yy and vertical velocity for June and December from 
the AER model. The parame erized K yy values shown are constant for six months of the 
year, changing in April and October. 

The LLNL model shows the strongest Hadley circulation in the troposphere. The GSFC 
model has a weak Hadley circulation, and the AER model has none. The GSFC circu- 
lation is the weakest overall n the stratosphere. The AER model has consistently the 
largest diffusion coefficients. Mean age of air calculated from the models reflects both 
advective and diffusive transport and can be compared with age fields derived from ob- 
servations. The GSFC model has the longest age of air in the upper stratosphere, with 
LLNL significantly shorter and AER shortest as reported in the M&M II report [ Park et 
ai, 1999]. Since most 2-D models generate age of air which is shorter than observations 
indicate, the GSFC model is most realistic in this regard. 

2.2. Diurnal Treatment and Time Stepping Scheme 

The approaches for dealing with time stepping and diurnal variability differ in significant 
ways among the three model. 1 . The AER and GSFC models use a family approach and 
transport only long-lived spedes, whereas the LLNL model does explicit time stepping 
of all species. The GSFC model transports the species N 2 Or> and C10N0 2 which the 
AER model considers part of the N0 y and/or Cl y families. Since zonal mean quantities 
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are transported in ‘2-D models, the zonal mean production and loss rates are used in the 

mass-continuity equation for long-lived species. Different techniques are used to obtain 
the zonal-mean production i nd loss rates. These involve first calculating the diurnal 
variations in the radical species. 

Computation of radical sp< cies in the AER model is performed with 17 explicit time 
steps over the diurnal cycle. Production and loss rates of long-lived species are com- 
puted from the radical concentrations at each of the 17 time steps and averaged over 24 
hours. New production and loss rates are computed daily. The Smolarkiewicz scheme 
[Smolarkiewicz, 1984] is used for time stepping of advection. The rates of change due to 
chemistry, diffusion, and advection are used to update long-lived species with a time step 
of 6 hours using an explicit scheme. 

The GSFC model does explicit diurnal calculations only for those radical species which 
are produced at night, ie. 1 NO3, NO2, NO3, N2O5, HOC1, HC1, CIO, CIONO2, BrO, 
HOBr, and BrON0 2 . For other species, daytime average photolysis rates are used, along 
with dawn values of the night species, to calculate daytime average radical concentrations 
and production and loss rate; for long-lived species. Time stepping of long-lived species 
is done with split operators, using a 12 hour time step for advection, a 3 hour step for 
vertical diffusion, a 24 hour si ep for horizontal diffusion, and a 24 hour step for chemistry. 
The Lin and Rood scheme [ L in and Rood, 1996] is used for advective transport. 

The LLNL model does exj -licit diurnal chemistry for all species with 15 minute time 
steps for chemistry and 2 hours time steps for transport, time marching 2 days at a time 
for each process. The Smolarkiewicz scheme is used for advection. To save computer time, 
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the diurnal version is run for one year while saving diurnal averaging coefficients daily. 

Then the diurnal average vers on of the model is run for two years using these coefficients. 

2.3. Chemistry Content 

Each model uses the JPL-97 compendium [DeMore et al. . 1997] as the source of reaction 
rate data. The monthly mean zonal mean temperature from climatology is used to cal- 
culate the rate constants for $as-phase reactions. Six heterogeneous reactions on sulfate 
aerosol are employed in each Model. Aerosol surface area is specified in these calculations, 
so reaction rates on sulfate aerosol are the same for all three models except for the way 
temperature variation is included in the rate calculation. Both LLNL and AER utilize 
a temperature distribution taxen from NCEP reanalysis data [Kalnay et al., 1996] which 
accounts for longitudinal and day-by-day deviations of temperature from the monthly 
zonal mean temperature. The reaction probability and the molecular thermal velocity are 
evaluated at each temperatun • in the distribution and then the mean value of the product 
is used, along with the aerosol surface area, to obtain the zonal mean reaction rate. This 
methodology is found to mak< • a significant difference in reaction rates for those reactions 
which are strongly temperature dependent (CIONO 2 + H 2 O, Br0N02 + H 2 O, CIONO 2 
+ HC1) [ Weisenstein et al., 1998; Pitari, 1993]. The GSFC model employs only the zonal 
mean temperature in the rate calculation for sulfate aerosols. 

Various intercomparison exercises have demonstrated that standard application of 
Beer’s law to account for 0 2 and 0 3 absorption in photolysis calculations provides re- 
liable results in most cases. ( 'are must be used, however, in the Schumann-Runge (S-R) 
bands for 0 2 photolysis and MO absorption. The AER model uses the photolysis code of 
Michael Prather [ Prather , 1903] which handles the S-R bands according to Minschwaner 
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et al. [1992], The online phot olysis calculation includes the effects of Rayleigh scattering 
in a spherical atmosphere. T1 e attenuated flux in 77 wavelength bands at each grid point 
is calculated and multiplied by the molecular cross sections obtained from DeMore et al. 
[1997]. Photolysis rates or J-rutes are the integral of this product over wavelength, and are 
calculated for the ten daytime time points. Photolysis rates for O 2 and NO are calculated 
separately to account for the ine structure of the S-R bands. 

The GSFC model uses a lookup table for the photolytic source term (PST) and pho- 
tolysis of O 2 [J(02)]. This table was generated by R. Kawa (GSFC) using a radiation 
code developed by D. Anderson and coworkers at the Johns Hopkins University Applied 
Physics Laboratory [Anderso 1 and Meier, 1979; Anderson and Lloyd, 1990]. The PST 
and J (O 2 ) table is given on tl e pressure grid of the GSFC 2-D model and is a function of 
wavelength, solar zenith angle and column ozone. The PST is computed for the particular 
wavelength, solar zenith angl \ column ozone, and pressure of interest and is multiplied 
by the solar irradiance at the 1 op of the model to compute the flux at the point of interest. 
The flux is then used with th< photodissociation cross section for a particular constituent 
to derive the photolysis rate < >f interest. 

The lookup table used in t he LLNL model is a derivative of the GSFC lookup table 
described above. The PST was generated at each of 79 wavelengths for LLNL by R. 
Kawa (GSFC) using the radiative transfer code mentioned above for 35 pressure levels, 
20 solar zenith angles, and 12 overhead ozone columns with varying range appropriate to 
the associated pressure. Molecular cross sections for photolysis for the same wavelength 
bin structure were determine d from laboratory sources [DeMore et al., 1997], at 200 
values of local temperature, i-overing the range encountered in the atmosphere. The J 
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values are then determined on-the-fly by integrating, over wavelength, the product of the 
exoatmospheric flux, the interpolated PST, and the actinic cross section at the current 
local temperature. Because evaluation of . 1 ( 0 - 2 ) requires treatment of the S-R bands, 
LLNL interpolated J(0 2 ) from the GSFC values as a function of pressure, solar zenith 
angle, and overhead ozone column, rather than by integration over the wavelength bins 
in the model. 

All three models include lightning as a source of N0„ in the tropical troposphere, with 
the AER and GSFC models as suming a source strength of 2 megatons per year distributed 
from 4 to 14 km in altitude. The LLNL model assumes a source strength of 5 megatons 
per year distributed according to the ISCCP cloud database [Rossow and Schiffcr, 1999]. 
Some species are removed in lie troposphere through rainout/ washout processes in each 
of the models, though both the rates of these processes arid the species subject to washout 
vary somewhat from model t > model. All three models remove H 2 0 2 , HN0 3 , HC1, and 
HBr. In addition to the common set, the GSFC model also removes H0 2 N0 2 , CH 3 OOH, 
HF, CC1F0, and CF 2 0; the AER model also removes CH 3 OOH, CH 3 O 3 H, CH 2 0, HF, and 
CF 2 0; and the LLNL model ilso removes N 2 0 5 , H0 2 N0 2 , CH 2 0, CIO, H0C1, C10N0 2 , 
and Br0N0 2 . Washout rates in the LLNL model vary from 3 days in the lowest 4 km to 
50 days near the tropopause The AER model uses similar washout rates (5 days near 
the surface to 40 days at 10 cm), but with no washout above 10 km. The GSFC model 
used the slowest washout rates: from 25 to 100 days. 

Distributions of H 2 0 in the troposphere are based on prescribed relative humidity values 
in the AER and LLNL models. In addition, the AER model specifies the water vapor 
mixing ratio just above the tiopopause (2.75 ppmv in the tropics, up to 3.5 ppmv at high 
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latitudes). The GSFC model ipecifies water vapor below 400 mb based on the climatology 
of Oort [1983]. Water vapor is transported in the troposphere above 400 mb in the GSFC 
model, with H 2 0 concentrati ms in excess of 50% relative humidity rained out [Fleming 
et al . , 1995]. All three mocels calculate H 2 0 concentration in the stratosphere using 
appropriate chemical sources {including methane oxidation) and sinks. 

2.4. Polar Heterogeneous Chemistry Formulation 

Though gas-phase reaction rates of stratospheric interest have been fairly well standard- 
ized by the JPL compendium DeMore et al., 1997], treatment of heterogeneous chem- 
istry, particularly under cold polar conditions, differs greatly between models. Approaches 
include assuming nitric acid rihydrate (NAT), sulfuric acid trihvdrate (SAT), or super- 
cooled ternary solution (STS) compositions for type I polar stratospheric cloud (PSC) 
particles. Some models calculate PSC surface area based on model calculated HN0 3 
and H 2 0, while others base PSC surface area on observations. Some models account for 
sedimentation of PSC particles, while others assume a fixed rate of denitrification under 
certain conditions. Two-dimensional models in addition need to account for the zonal 
asymmetry of temperature. Since there is no consensus on how best to represent cold 
polar processes, the GSFC, I LNL, and AER models implement different approaches. 

Both GSFC and AER employ thermodynamic equilibrium parameterizations to calcu- 
late Type I (assumed to be NAT) and Type II (ice) PSCs based on available gas phase 
H 2 0 and HNO3 and a distribution of temperatures. The GSFC model assumes that su- 
persaturation factors of 10 for NAT and 1.4 for ice are required before PSCs form. The 
AER model assumes no supe: saturation. The GSFC model assumes lognormal size distri- 
butions of PSC particles, with NAT having a mode radius of 1.0 fim and a a of 1.8 and ice 
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particles having a mode radius of 10 /mi and a a of 1.8 [Considine et al., 1994], The AER 

model assumes a single radii s for each type of particle. 0.5 pm for NAT and 7 pin for 
ice. Because of the particle si/.e assumptions, the GSFC model has greater sedimentation 
and denitrification/dehydration. The temperature distribution is applied only to obtain 
the PSC surface area in the ( ISFC model, while the AER model applies the temperature 
distribution to obtain the pr< duct of the reaction rate and surface area. 

The LLNL model handles cold heterogeneous processes by assigning a supercooled 
ternary solution (STS) composition to PSC particles. The particle surface area is not 
calculated by the model. Instead, a surface area density of 1 pm 2 /cm 3 is imposed within 
25° of the poles when the PSC climatology of Poole and Pitts [1994] indicates a PSC 
frequency of occurrence exceeding 0.08. Reaction rates are obtained by integrating over 
the temperature distribution and then multiplying by the fixed PSC surface area. Dehy- 
dration and denitrification ar ; represented globally (independent of the PSC surface area 
parameterization) by assuming that the partial pressure in excess of the saturation va- 
por pressure over ice (calcula' ed using zonal mean temperatures) is removed permanently 
with first order time constants of 1 day for H 2 O and 0.5 days for HNO 3 . 

3. The Hybrid Models 
3.1. Approach 

Transport parameters from the GSFC and LLNL models were used in the AER model. 
The version of the AER model with GSFC transport coefficients will be referred to as 
the AER/GSFC model. The version of the AER model with LLNL transport coefficients 
will be referred to as the AER/LLNL model. The streamfunction can be interpolated to 
different grid resolutions whi e maintaining a nondivergcnt flow field. The horizontal and 
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vertical diffusion fields, K yy and K 22 , are also used and can be interpolated between grids 
if necessary. The GSFC and LLNL models do not use a K yz field, so I< y2 is set to zero in 
the hybrid models. 

The transport algorithms and numerics of the AER model were modified to some extent 
for the hybrid models. It w is found necessary to use the LLNL grid (5° latitude, 1.5 
km vertical) within the AER /LLNL model when using the LLNL transport parameters 
because of sharp gradients i i the LLNL diffusion fields. The AER/LLNL model was 
also modified to use the sam ■ time splitting scheme for the advective and diffusive time 
steps as the LLNL model. Since the AER and LLNL models both use the Smolarkiewicz 
transport scheme, the transport algorithm was not modified. The top boundary of the 
AER/LLNL model was set a’ 80 km as in the LLNL model. 

The GSFC model uses the Lin and Rood [1996] scheme for transport. We decided that 
it was not practical to modib the transport scheme for the AER/ GSFC model, so it uses 
Smolarkiewicz like the AER model. We ran the AER/GSFC model in the original GSFC 
grid resolution and in the AE R grid resolution. In both cases, differences with the GSFC 
model were found due to tb' different transport scheme, and the grid resolution made 
only a small additional difference. Therefore we opted to employ the AER model grid 
resolution in the AER/GSFC results presented here. Note that the top boundary of the 
AER/GSFC model is 60 km. not 90 km as in the GSFC model. 

The hybrid models use the chemistry formulation of the AER model, including the 
photolysis scheme, reaction r; tes, and heterogeneous chemistry formulation for sulfate and 
PSCs. The family approach and diurnal calculations of the AER model are retained in 
the hybrids. The hybrid models thus represent the AER model s chemistry with dynamics 
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replaced by those of another nodel. The following inert tracer experiment is used to test 
whether the LLNL and GSFO transport formulations have been implemented correctly 
within the hybrid models, and how closely the implementation matches the native model. 

3.2. Inert Tracer Comparison 

The Models and Measurements II experiment A-3 [ Park et al., 1999] is used to test the 
accuracy of the GSFC and LLNL transport implementation within the AER model. This 
experiment is designed to mimic the NO y perturbation due to a fleet of 500 supersonic 
aircraft cruising at Mach 2.4. with most emissions occurring in the northern hemisphere 
at altitudes of 18-20 km. The source is taken from the 1995 NASA HSCT scenario 
[Stolarski et al, 1995] with NO x emissions of 10 grams N0 2 per kilogram of aircraft fuel 
burned. Removal in the troposphere is simulated by setting a boundary condition of zero 
concentraiton below 6 km. T lere is no stratospheric removal. 

June results of the inert tracer calculation from the AER, AER/GSFC, AER/LLNL, 
GSFC, and LLNL models are shown in Figure 4. Comparison of the results from the 
native models (panels a, d, and e) shows that the GSFC model retains the most emitted 
material in the stratosphere (80% more than the AER model and 65% more than the 
LLNL model). Consequently the GSFC model transports the most tracer to the upper 
stratosphere and the southern hemisphere (by more than a factor of two). When these 
results are scaled to represent equal atmospheric burdens of injected tracer, the GSFC 
model obtains 18% more tra< er at 45°S and 20 km than the AER model and 29% more 
than the LLNL model, indie iting that, the models differ in their global distribution of 
tracer. Differences between the native model and the corresponding hybrid model reflect 
the extent to which the hybrid model transport does not reproduce the native model. 
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Differences between the GSF( ' and AER/GSFC models are 4-6% above 20 km and much 
larger near the tropopause (ip to a factor of two). The comparison between the LLNL 
and AER/LLNL models shov s differences of only 1-3% everywhere. 

Differences between the GSFC and AER/GSFC models are likely due largely to the 
difference in the transport schemes, Smolarkiewicz [1984] for the AER model and Lin 
and Rood [1996] for the GSFG model. The difference in transport schemes is especially 
noticeable near the tropopai se where gradients are large. Note that the AER/GSFC 
model results were obtained with the GSFC transport parameters interpolated to the 
AER grid. Performing the cal mlation on the GSFC grid in the AER/GSFC model results 
in only slightly better agreement. The LLNL and AER/LLNL models both use the Smo- 
larkiewicz transport scheme and show minor differences. We attribute these differences 
to minor differences in the model numerics and the Smolarkiewicz implementation. This 
set of exercises assures us that the hybrid models provide a good approximation of the 
transport properties of the native models. 

4. Background Atmosphere Intercomparison 

We now use the results from the hybrid models to illustrate how they can help to explain 
the differences in simulated distributions of chemical species in the atmosphere. The 
background atmosphere is obtained from the Atmospheric Effects of Aviation Program 
(AEAP) calculations that coriesponds to a 2015 atmosphere with subsonic aircraft [ IPCC , 
1999; Kawa et al . , 1999], and is labeled scenario D in the IPCC [1999] calculations. The 
2015 atmosphere is assumed to contain a stratospheric Cl,, concentration of 3.0 ppbv, N 2 0 
surface concentration of 330 >pbv, and background (nonvolcanic) levels of stratospheric 


aerosol. 
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Each of the calculations pi esented represents the model’s steady-state condition, i.e. 
there is no appreciable change in the species’ concentrations from one year to the next. 
Because of the uncertainties associated with PSC chemistry, two simulations were per- 
formed for each model: one \\ ith and one without PSC chemistry. The model results are 
discussed in the following foimat. Figure 5 shows the distributions of NO y , H 2 0, Cl y , 
and 0 3 from the AER model in October with the PSC parameterization. The differences 
among models are discussed using Figures 6-15. In these figures, panels (a), (b) and (c) 
show the differences in calcul ited concentrations for cases with a PSC parameterization 
between the models for (X-AER), (AER/X-AER), and (X-AER/X) respectively, where X 
represents either GSFC or LLXL and AER/X is a hybrid model. Note that algebraically, 
the sum of the differences in panels (b) and (c) equals the difference in panel (a). The 
difference in panel (b) is due ;.o differences in transport parameters (winds and eddy dif- 
fusion coefficients). The difference in panel (c) is mostly due to differences in chemical 
formulation. However, differences in numerical treatment and time-stepping also con- 
tribute. Panel (d) shows the difference (X-AER/X) without the PSC parameterization. 
Comparison of panels (c) and (d) isolates the difference due to PSC treatment. 

It is not the purpose of this-, study to compare the model results with observations. In 
fact, use of the 2015 background atmosphere makes observational comparisons impossible. 
It is useful to keep in mind, however, that the AER model tends to calculate too low a 
concentration in the lower stratosphere for downward diffusing species (eg. 0 3 , NO y , Cl y ). 
Compared with the TOMS climatology, the AER model calculates too much ozone during 
the springtime maximum (try about 10%-15%), and underpredicts ozone column at the 
tropics by about 5%-10%. 
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4.1. NO y 

Figure 5a shows the distribution of total odd nitrogen (NO y ) calculated by the AER 
model for October. A maximum value of 20 ppbv is calculated in the tropics at around 
40 km, falling off to 10 ppb\ at 50kin and 30 km. At higher latitudes, the maximum 
NO y value (about 15 ppbv) is found near 30 km. The GSFC and LLNL native models 
(Figures 6a and 7a) calculate more NO y than the AER model (by as much as 12 ppbv 
in the upper stratosphere). Transport differences (shown in Figures 6b and 7b) lead to 
relatively large differences (both positive and negative, up to 6 ppbv) in concentration in 
the upper stratosphere and tiopical mid stratosphere as well. With the premise that the 
hybrid models mostly remov< the effects of different transport, the differences shown in 
Figures 6c and 7c should be ihe result of different chemistry. For NO y , this would mean 
the source terms from N 2 0 and lightning, the gas-phase removal term by photolysis of 
NO, removal by PSC chemisti y, and washout in the troposphere. Results from Figures 6d 
and 7d show that the effect oi PSCs is limited to near the poles at 10-25 km altitude. The 
upper stratospheric chemical difference is likely due to differences in NO photolysis. 

4.2. H 2 0 

Figure 5b shows the gas-phase H 2 0 distribution calculated by the AER model for Oc- 
tober. Figure 8a shows that i he GSFC model calculates up to 1 ppmv more H 2 0 in the 
lower and middle stratosphere. Figures 8b and 8c suggest that most of the difference in 
the lower stratosphere is due to the difference in the tropospheric boundary condition. 
However, the transport differ* nee is still responsible for half the difference in some regions 
of the middle stratosphere. T ie LLNL model also calculates more H 2 0 (Figure 9a) in the 
lower stratosphere and much of the middle stratosphere. In this case, transport differences 
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are generally larger than the chemical differences, but of opposite sign. Figures 8d and 

9 d show that differences in PSC treatments have a small effects in both cases. 

4.3. Cl y 

The distribution of total inorganic chlorine (Cl ;/ ) calculated by the AER model for Oc- 
tober is shown in Figure 5c. The differences among the native models (Figures 10a and 
11a) in the upper stratosphere are less than 5%. This shows that each model achieves an 
inorganic chlorine content in i he upper stratosphere roughly equal to the input of organic 
chlorine in the troposphere. 1 he differences are much larger in the lower stratosphere, ap- 
proaching a factor of 2. Transport is responsible for the bulk of the differences (Figures 10b 
and lib). However, Figures LOc and 11c still show differences due to chemistry, which 
can be attributed to differem es in washout and details of the numerical schemes (which 
affect the simulation of transport in regions with large gradients such as the tropopause). 
PSC treatments have a negligible effect on Cl y . 

4.4. Ozone 

Figure 5d shows the O3 distribution calculated by the AER model for October. The 
native GSFC (Figure 12a) and LLNL (Figure 13a) models show less 0 3 (by as much as 
15% or 1.5 ppmv) than the AER model in the equatorial middle stratosphere and the 
midlatitude lower to middle stratosphere. About 2/3 of the difference in the tropics is 
due to transport and the renninder due to chemistry (see Figures 12b, c and 13b, c). In the 
southern high latitudes in Oc ober, chemistry accounts for most of the difference between 
the AER and GSFC models, s pecifically the PSC treatment (as can be seen by comparing 
Figures 12c and d). The comparison between the AER and LLNL models in the southern 
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high latitudes shows that effects of transport and chemistry partially cancel each other in 
October, with transport having a somewhat larger effect. 

The results from comparisons of calculated ozone column are given in Figures 14 and 
15. Both the GSFC and the I LNL models calculate more 0 3 in the tropics and less 0 3 at 
high latitudes than the AER model (generally resulting in better agreement with TOMS 
observations). The LLNL model also calculates more 0 3 in the midlatitudes than the AER 
model. Both chemistry and t ransport play significant roles in the differences, often with 
competing effects. In the soul hern hemisphere high latitudes, the GSFC model calculates 
20-100 DU less than the AER model. Transport differences between the AER and GSFC 
models produce more ozone in austral summer and less in austral winter. In the absence 
of PSCs (see Figure 14d), chemical differences account for a difference of less than 20 DU. 
Inclusion of PSCs in the models greatly reduces the ozone column in the GSFC model 
relative to the AER model. In the AER-LLNL model comparisons, transport leads to 
lower ozone columns at high southern latitudes for LLNL, while chemistry leads to higher 
ozone columns in the same region. Without PSCs in either model, chemical differences 
are less than 10 DU. The LLNL model has the least efficient PSC mechanism for depleting 
ozone and the GSFC model the most efficient scheme. 

5. Intercomparison of Calculated HSCT Perturbations 

Over the past few decades there has been much interest in a potential fleet of commercial 
supersonic aircraft and its atmospheric impact. There have been numerous modeling 
studies to investigate the impact of High Speed Civil Transport (HSCT) aircraft in the 
lower stratosphere and upper troposphere [ IPCC , 1999; Kawa et al., 1999; and references 
therein]. These studies consist of various scenarios used to explore the effects of parametric 
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modifications in the HSCT aircraft fleet size, cruise altitude, emission parameters, and 
the background atmosphere. The model response can be generally explained in terms 
of the amount of engine emissions retained in the stratosphere (a transport issue), and 
the response of ozone to thes< perturbations (where both transport and chemistry play a 
role). It is the latter complic ition that makes it difficult to separate the two effects and 
intercompare model results in a meaningful way. 

To make use of existing computations, we used results from the IPCC [1999] exercise, 
in which the AER, GSFC and the LLNL models participated, and generated additional 
results with the AER/GSFC and the AER/LLNL hybrid models. We have chosen to 
examine scenario Sic from iPCC [1999], which includes a fleet of 500 HSCT aircraft 
operating in 2015, with a N(). r emission index of 5 grams of NO 2 per kilogram of fuel 
burned and H 2 O emission index of 1230 grams per kilogram of fuel. The aircraft have 
a cruise altitude of 18-20 km and fly predominately in the northern hemisphere. This 
scenario is compared with Scenario D discussed in section 4 which represents a 2015 
background atmosphere with only subsonic aircraft. 

Scenario Sic was calculated by nine modeling groups for the IPCC [1999] assessment 
report. Annual average ozone column perturbations ranged from 0.0% to -0.4%' for the 
northern hemisphere and from 0.0% to -0.8% for the southern hemisphere [IPCC, 1999]. 
This large difference in ozone response is difficult to interpret without a complete under- 
standing of the model differences. By examining perturbations of 0 3 column, local 0 3 , 
H 2 0, and NO y using the AER/GSFC and the AER/LLNL hybrid models, we can ob- 
tain an understanding of how the different model chemical formulations affect the ozone 
perturbations. This, in turn may provide a better understanding of the variability in 
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the calculated ozone perturbations from different modeling groups. We will first discuss 
calculations without PSC chemistry, to simplify the interpretation. We will then show 
O 3 perturbation results with PSC chemistry, which has quite a large impact on these 
calculations. 

Perturbations to H 2 O are shown in Figure 16 for October without PSC chemistry. 
Perturbations to NO y are shown in Figure 17 for the same conditions. The H 2 0 results 
are similar to those for the inert tracer shown in Figure 4 when scaled for emissions (factor 
of 0.31). The N() y results should scale by 0.5 from results shown in Figure 4 except that 
NOy has chemical loss in th * upper stratosphere. As in experiment A-3, more of the 
emissions are transported to he southern hemisphere in the GSFC model than the AER 
or LLNL models. The HSCT NOy and 11 2 0 perturbation differences between the native 
models and the hybrid model; are fairly small, despite the fairly large differences between 
the NOy and H 2 0 backgroun 1 concentrations, indicating that transport plays the major 
role for these gases. 

Perturbations in O 3 for 0< tober due to HSCT emissions are shown in Figure 18 for 
simulations without PSC chemistry. The GSFC and AER/GSFC models show twice as 
much O 3 depletion in the upper stratosphere as the other models due to having twice the 
H 2 0 perturbation in this region. Models with similar transport show similar ozone deple- 
tion in the upper stratosphere, but somewhat different ozone perturbations in the lower 
stratosphere. These different'' >s could be due to the different background levels of NOy, 
Cly, Br y , and H 2 0 in this region or to different chemical treatment of 0 3 . Differences in 
transport (compare panels a, b, and c) have a large impact on the calculated O 3 pertur- 
bation near the tropopause al mid and high latitudes. Figure 19 shows perturbations in 
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0 3 column for simulations wit bout PSC chemistry. Of the native models, AER and LLNL 
show only 0.1% O3 depletion in the southern tropics, while the GSFC model shows 0.3% 
O3 depletion here. The LLNL model shows the greatest 0 3 depletion at both poles. The 
GSFC model shows the least depletion at the north pole, but the most at 30°N. 

The seasonal variability at 1 igh latitudes seen in the AER model is reflected in the hybrid 
models as well, indicating that this feature is due to the chemical scheme in the AER 
model and not due to transport. The LLNL model calculates much more O3 depletion in 
the high southern latitudes t lan the AER/LLNL model. Background concentrations of 
NOy in the high latitude lower stratosphere are much larger in the LLNL model than the 
AER/LLNL model and background H 2 0 concentrations are lower, which would result in 
more sensitivity to NO,, perturbations in the LLNL model. The GSFC and AER/GSFC 
models both show a maxima of O3 depletion at 30°N, indicating that this is a feature of 
the GSFC transport (a result of the O3 enhancement from 40-90°N at 8-12 km altitude 
seen in Figure 18). However, the AER/GSFC model shows a maximum of O3 depletion 
at the north pole, especially in springtime, whereas the GSFC model shows a minimum 
here. 

With PSC chemistry, all me <dels show denitrification and dehydration at the south pole 
in austral spring. The GSFC model, but not the AER or LLNL models, also shows 
denitrification near the north pole in springtime. This is due to the different ways that 
NAT particle sizes are specified in the models, resulting in more particle sedimentation 
in the GSFC model. O3 perturbations for October with PSC chemistry are shown in 
Figure 20, and perturbations in O3 column are shown in Figure 21. All models show 
enhanced O3 depletion in both hemispheres with PSC chemistry included, except for the 
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LLNL model in the northeri hemisphere. PSC chemistry has a very large impact on 
calculated ozone perturbations due to HSCT at high latitudes and increases the model 
differences. The hybrid models show much more 0 3 depletion at high northern latitudes 
than the GSFC and LLNL models, though the GSFC model shows more depletion than 
the AER/GSFC model in the southern high latitudes. 

6. Discussion and Conclusions 

Traditional comparisons oi model-calculated ozone perturbations are incomplete be- 
cause reasons for the model differences cannot be diagnosed. We have approached this 
problem by constructing hybrid models which employ the transport fields of one model 
within the framework of another model. We demonstrate that the use of hybrid models 
can help to separate, to a large extent, the effects of transport and chemistry on model cal- 
culated results. The intercom parison performed in this study has lead to improvements in 
the participating models. Because of the NO ?/ comparison problem identified in this study, 
the LLNL model has modified J(NO) by creating a new lookup table for this parameter 
derived from AER results (ie. the UCI photolysis code) calculated for January and July, 
used by interpolation to the appropriate pressure, solar zenith angle, and overhead ozone 
column. 

The inert tracer results in section 3.2 show that the native models’ advective winds 
and eddy coefficients accoun s for the bulk of the transport differences, although the 
numerical scheme still has a small effect. The comparisons of the background atmosphere 
indicate that the tropospheric boundary condition for H 2 0 has a large effect on the model 
calculated stratospheric distribution. The very large difference in NO y in the middle and 
upper stratosphere can probably be attributed to treatment of the NO photolysis rate. 
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The choice of transport parameters has a substantial effect on model calculated NO y , Cl y 
and H 2 0 in the lower and middle stratosphere, with both positive and negative differences 
seen. Differences in ozone an the result of differences in NO,,, Cl y , Br y , and H 2 0 caused 
by both chemistry and transport, as well as differences in transport of ozone itself. PSC 
parameterizations represent < large difference in model formulations, and this is evident 
in the ozone comparisons. 

In the HSCT perturbation calculations, our results showed that most of the differences 
in the H 2 0 and NO y perturbations are due to transport. Calculated ozone differences in 
the upper stratosphere can be attributed to differences in transport of H 2 0. Differences 
in the lower stratosphere are nore difficult to interpret, but high latitude ozone behavior 
appears to be controlled by chemistry to a greater extent than midlatitude ozone. When 
PSC processes are omitted from the models, differences in ozone column perturbation 
between models with the sam ' transport are fairly small (less than 0.2% in most regions). 
When PSC processes are included, differences can be as high as 1% at high latitudes. 
For models with identical chemistry but different transport, differences in ozone column 
perturbation due to HSCT ere no larger than 0.3% when PSC processes are omitted. 
Including PSCs yields differences in ozone column perturbation of up to 1.6% between 
models with identical chemist ry. 

This study has allowed us t< quantitatively diagnose transport differences between three 
of the models used in the II CC assessment of supersonic aircraft [IPCC. 1999], It is 
well known that differences in model transport play a major role in calculated ozone 
perturbation differences. A more surprising result is that differences in model chemical 
treatment produce significant differences in perturbation results as well, despite the fact 
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that all models use the same set of reaction rate coefficients. These differences are much 

larger when PSC processes are included in the models. As more is learned about the 
formation of PSCs in the futvre. it is hoped that model treatment of these processes will 
become more standardized and this source of model discrepancy reduced. 
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Figure 1 . LLNL model val les of K yy (in 10 10 cm 2 /s) (panels a and b) and vertical velocity (in 
mm/s) (panels c and d) for June (panels a and c) and December (panels b and d). 

Figure 2. GSFC model values of K yy (in 10 10 cm 2 /s) (panels a and b) and vertical velocity 
(in mm/s) (panels c and d) for June (panels a and c) and December (panels b and d). 

Figure 3. AER model values of K yy (in 10 10 cm 2 /s) (panels a and b) and vertical velocity (in 
mm/s) (panels c and d) for June (panels a and c) and December (panels b and d). 

Figure 4. Calculated mix ng ratio (ppbv) for June of an inert tracer with source similar to 
HSCT aircraft emission of N( > y and tropospheric sink. Panel (a) shows AER model results, panel 
(b) AER/GSFC results, panel (c) AER/LLNL results, panel (d) GSFC results, and panel (e) 
LLNL results. 

Figure 5. October mixing ratios of (a) N0 y in ppbv, (b) H 2 O in ppmv, (c) Cl y in ppbv, (d) 
O3 in ppmv, and (e) O 3 column in dobson units as calculated by the AER model for steady-state 
2015 conditions. 

Figure 6. October NO y in xing ratio differences (ppbv) between the GSFC and AER models 
in a 2015 background atmosphere with subsonic aircraft. Panel (a) shows the difference between 
the GSFC and AER models with PSC chemistry included, panel (b) the difference between the 
AER/GSFC and AER models, and panel (c) the difference between the GSFC and AER/GSFC 
models. Panel (d) shows the difference between the GSFC and AER/GSFC models without PSC 


chemistry. 
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Figure 7. October N() y mixing ratio differences (ppbv) between the LLNL and AER models 
in a 2015 background atmosphere with subsonic aircraft. Panel (a) shows the difference between 
the LLNL and AER mode s with PSC chemistry included, panel (b) the difference between the 
AER/LLNL and AER models, and panel (c) the difference between the LLNL and AER/LLNL 
models. Panel (d) shows the difference between the LLNL and AER/LLNL models without PSC 
chemistry. 


Figure 8. October H 2 0 mixing ratio differences (ppmv) between the GSFC and AER models 
in a 2015 background atmosphere with subsonic aircraft. Panel (a) shows the difference between 
the GSFC and AER models with PSC chemistry included, panel (b) the difference between the 
AER/GSFC and AER models, and panel (c) the difference between the GSFC and AER/GSFC 
models. Panel (d) shows the difference between the GSFC and AER/GSFC models without PSC 
chemistry. 


Figure 9. October H 2 0 mixing ratio differences (ppmv) between the LLNL and AER models 
in a 2015 background atmosphere with subsonic aircraft. Panel (a) shows the difference between 
the LLNL and AER models with PSC chemistry included, panel (b) the difference between the 
AER/LLNL and AER models, and panel (c) the difference between the LLNL and AER/LLNL 
models. Panel (d) shows the difference between the LLNL and AER/LLNL models without PSC 


chemistry. 
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Figure 10. October Cl y mixing ratio differences (ppbv) between the GSFC and AER models 
in a 2015 background atmosphere with subsonic aircraft. Panel (a) shows the difference between 
the GSFC and AER models with PSC chemistry included, panel (b) the difference between the 
AER/GSFC and AER models, and panel (c) the difference between the GSFC and AER/GSFC 
models. Panel (d) shows the difference between the GSFC and AER/GSFC models without PSC 
chemistry. 


Figure 11. October Cl y mixing ratio differences (ppbv) between the LLNL and AER models 
in a 2015 background atmosphere with subsonic aircraft. Panel (a) shows the difference between 
the LLNL and AER models -with PSC chemistry included, panel (b) the difference between the 
AER/LLNL and AER models, and panel (c) the difference between the LLNL and AER/LLNL 
models. Panel (d) shows the difference between the LLNL and AER/LLNL models without PSC 
chemistry. 


Figure 12. October O 3 mixing ratio differences (ppmv) between the GSFC and AER models 
in a 2015 background atmosphere with subsonic aircraft. Panel (a) shows the difference between 
the GSFC and AER models -with PSC chemistry included, panel (b) the difference between the 
AER/GSFC and AER models, and panel (c) the difference between the GSFC and AER/GSFC 
models. Panel (d) shows the difference between the GSFC and AER/GSFC models without PSC 


chemistry. 
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Figure 13. October O 3 mixing ratio differences (ppmv) between the LLNL and AER models 
in a 2015 background atmosphere with subsonic aircraft. Panel (a) shows the difference between 
the LLNL and AER models with PSC chemistry included, panel (b) the difference between the 
AER/LLNL and AER models, and panel (c) the difference between the LLNL and AER/LLNL 
models. Panel (d) shows ti e difference between the LLNL and AER/LLNL models without PSC 
chemistry. 

Figure 14. O 3 column differences (Dobson units) between the GSFC and AER models in 
a 2015 background atmosphere with subsonic aircraft. Panel (a) shows the difference between 
the GSFC and AER mode s with PSC chemistry included, panel (b) the difference between the 
AER/GSFC and AER models, and panel (c) the difference between the GSFC and AER/GSFC 
models. Panel (d) shows the difference between the GSFC and AER/GSFC models without PSC 
chemistry. 

Figure 15. O 3 column differences (Dobson units) between the LLNL and AER models in 
a 2015 background atmosphere with subsonic aircraft. Panel (a) shows the difference between 
the LLNL and AER mode's with PSC chemistry included, panel (b) the difference between the 
AER/LLNL and AER models, and panel (c) the difference between the LLNL and AER/LLNL 
models. Panel (d) shows the difference between the LLNL and AER/LLNL models without PSC 
chemistry. 

Figure 16. Calculated perturbation in H 2 0 (ppmv) for October due to emission from HSCT 
aircraft in 2015 relative to an atmosphere with only subsonic aircraft. Panel (a) shows AER 
model results, panel (b) AER/GSFC results, panel (c) AER/LLNL results, panel (d) GSFC 
results, and panel (e) LLN L results, all without PSC chemistry. 
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Figure 17. Calculated perturbation in NO y (ppbv) for October due to emission from HSCT 
aircraft in 2015 relative to ail atmosphere with only subsonic aircraft. Panel (a) shows AER 
model results, panel (b) AER/GSFC results, panel (c) AER/LLNL results, panel (d) GSFC 
results, and panel (e) LLNL results, all without PSC chemistry. 

Figure 18. Calculated perturbation in 0 3 (%) for October due to emission from HSCT aircraft 
in 2015 relative to an atmosphere with only subsonic aircraft. Panel (a) shows AER model results, 
panel (b) AER/GSFC results, panel (c) AER/LLNL results, panel (d) GSFC results, and panel 
(e) LLNL results, all without PSC chemistry. 

Figure 19. Calculated perl urbation in 0 3 column (%) due to emission from HSCT aircraft in 
2015 relative to an atmosphere with only subsonic aircraft.. Panel (a) shows AER model results, 
panel (b) AER/GSFC results, panel (c) AER/LLNL results, panel (d) GSFC results, and panel 
(e) LLNL results, all without PSC chemistry. 

Figure 20. Calculated pert urbation in 0 3 (%) for October due to emission from HSCT aircraft 
in 2015 relative to an atmosphere with only subsonic aircraft. Panel (a) shows AER model results, 
panel (b) AER/GSFC results, panel (c) AER/LLNL results, panel (d) GSFC results, and panel 
(e) LLNL results, all with PSC chemistry. 

Figure 21. Calculated peri urbation in 0 3 column (%) due to emission from HSCT aircraft in 
2015 relative to an atmosphere with only subsonic aircraft. Panel (a) shows AER model results, 
panel (b) AER/GSFC results, panel (c) AER/LLNL results, panel (d) GSFC results, and panel 
(e) LLNL results, all with PSC chemistry. 





Figure 1 . LLNL model value; of K yy (in 10 10 cm 2 /s) (panels a and b) and vertical velocity (in mrn/s) 
(panels c and d) for June (panels a and c) and December (panels b and d). 
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Figure 2. GSFC model value; of K vv (in 10 10 cm 2 /s) (panels a and b) and vertical velocity (in mm/s) 
(panels c and d) for June (panels a and c) and December (panels b and d). 
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Figure 3. AER model values of K yy (in 10 10 cm 2 /s) (panels a and b) and vertical velocity (in mm/s) 
(panels c and d) for June (panels a and c) and December (panels b and d). 
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Figure 4. Calculated mixing ratio (ppbv) for June of an inert tracer with source similar to HSCT 
aircraft emission of N0 V and tropospheric sink. Panel (a) shows AER model results, panel (b) 
AER/GSFC results, panel (c) AER/LLNL results, panel (d) GSFC results, and panel (e) LLNL results. 
Contours are 0.01, 0.1, 0.2-2. 4 by 0.2, 2.8, and 3.2. 
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Figure 5. October mixing ra ios of (a) N0„ in ppbv, (b) H 2 0 in ppmv, (c) Cl„ in ppbv, (d) 0 3 in ppmv 
and (e) 0 3 column in dobson units as calculated by the AER model for steady-state 2015 conditions. 







Figure 6. October N0 y mixi ig ratio differences (ppbv) between the GSFC and AER models in a 2015 
background atmosphere with subsonic aircraft. Panel (a) shows the difference in the GSFC and AER 
models with PSC chemistry included, panel (b) is the AER/ GSFC- AER difference, and panel (c) shows 
the GSFC - AER/GSFC modi 1 difference. Panel (d) is the GSFC - AER/GSFC model difference without 
PSC chemistry. 
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Figure 7. October .\ ( ) v mixing ratio differences (ppbv) between the LLNL and AER models in a 2015 
background atmosphere with s-ubsonic aircraft. Panel (a) shows the difference in the LLNL and AER 
models with PSC chemistry in cluded, panel (b) is the AER/LLNL-AER difference, and panel (c) shows 
the LLNL - AER/LLNL model difference. Panel (d) is the LLNL - AER/LLNL model difference without 
PSC chemistry. 
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Figure 8. October H 2 O mixing ratio differences (ppmv) between the GSFC and AER models in a 2015 
background atmosphere with subsonic aircraft. Panel (a) shows the difference in the GSFC and AER 
models with PSC chemistry included, panel (b) is the AER/GSFC-AER difference, and panel (c) shows 
the GSFC - AER/GSFC model difference. Panel (d) is the GSFC - AER/GSFC model difference without 
PSC chemistry. 
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Figure 9. October H 2 O mixing ratio differences (ppmv) between the LLNL and AER models in a 2015 
background atmosphere with subsonic aircraft. Panel (a) shows the difference in the LLNL and AER 
models with PSC chemistry included, panel (b) is the AER/LLNL-AER difference, and panel (c) shows 
the LLNL - AER/LLNL model difference. Panel (d) is the LLNL - AER/LLNL model difference without 
PSC chemistry. 
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Figure 10. October Cl,, mixing ratio differences (ppbv) between the GSFC and AER models in a 2015 
background atmosphere with subsonic aircraft. Panel (a) shows the difference in the GSFC and AER 
models with PSC chemistry included, panel (b) is the AER/GSFC-AER difference, and panel (c) shows 
the GSFC - AER/GSFC mode I difference. Panel (d) is the GSFC - AER/GSFC model difference without 
PSC chemistry. 
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Figure 11. October Cl y mixing ratio differences (ppbv) between the LLNL and AER models in a 2015 
background atmosphere with subsonic aircraft. Panel (a) shows the difference in the LLNL and AER 
models with PSC chemistry included, panel (b) is the AER/LLNL-AER difference, and panel (c) shows 
the LLNL - AER/LLNL model difference. Panel (d) is the LLNL - AER/LLNL model difference without 
PSC chemistry. 
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Figure 12. October O 3 mixing ratio differences (ppmv) between the GSFC and AER models in a 2015 
background atmosphere with i-ubsonic aircraft. Panel (a) shows the difference in the GSFC and AER 
models with PSC chemistry included, panel (b) is the AER/GSFC- AER difference, and panel (c) shows 
the GSFC - AER/GSFC modt 1 difference. Panel (d) is the GSFC - AER/GSFC model difference without 
PSC chemistry. 
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Figure 13. October O 3 mixing ratio differences (ppmv) between the LLNL and AER models in a 2015 
background atmosphere with subsonic aircraft. Panel (a) shows the difference in the LLNL and AER 
models with PSC chemistry included, panel (b) is the AER/LLNL-AER difference, and panel (c) shows 
the LLNL - AER/LLNL model difference. Panel (d) is the LLNL - AER/LLNL model difference without 
PSC chemistry. 








Figure 14 . O 3 column differences (Dobson units) between the GSFC and AER models in a 2015 
background atmosphere with subsonic aircraft. Panel (a) shows the difference between the GSFC and 
AER models with PSC chemistry included, panel (b) the difference between the AER/GSFC and AER 
models, and panel (c) the difft rence between the GSFC and AER/GSFC models. Panel (d) shows the 
difference between the GSFC and AER/GSFC models without PSC chemistry. 
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Figure 15. O 3 column differences (Dobson units) between the LLNL and AER models in a 2015 
background atmosphere with subsonic aircraft. Panel (a) shows the difference between the LLNL and 
AER models with PSC chemistry included, panel (b) the difference between the AER/ LLNL and AER 
models, and panel (c) the difference between the LLNL and AER/ LLNL models. Panel (d) shows the 
difference between the LLNL and AER /LLNL models without PSC chemistry. 
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Figure 16 . Calculated perturbation in H 2 0 (ppmv) for October due to emission from HSCT aircraft in 
2015 relative to an atmosphere with only subsonic aircraft. Panel (a) shows AER model results, panel (b) 
AER/GSFC results, panel (c) AER/LLNL results, panel (d) GSFC results, and panel (e) LLNL results, 
all without PSC chemistry. 
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Figure 17. Calculated pertuibation in NO y (ppbv) for October due to emission from HSCT aircraft in 
2015 relative to an atmosphere with only subsonic aircraft. Panel (a) shows AER model results, panel (b) 
AER/GSFC results, panel (c) AER/LLNL results, panel (d) GSFC results, and panel (e) LLNL results, 
all without PSC chemistry. 
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Figure 18 . Calculated pertuibation in 0 3 (%) for October due to emission from HSCT aircraft in 2015 
relative to an atmosphere with only subsonic aircraft. Panel (a) shows AER model results, panel (b) 
AER/GSFC results, panel (c) AER /LLNL results, panel (d) GSFC results, and panel (e) LLNL results, 
all without. PSC chemistry. 
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Figure 19. Calculated perturbation in 0 3 column (%) due to emission from HSCT aircraft in 2015 
relative to an atmosphere with only subsonic aircraft. Panel (a) shows AER model results, panel (b) 
AER/GSFC results, panel (c) AER/LLNL results, panel (d) GSFC results, and panel (e) LLNL results, 
all without PSC chemistry. 
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Figure 20. Calculated perturbation in 0 3 (%) for October due to emission from HSCT aircraft in 2015 
relative to an atmosphere with only subsonic aircraft. Panel (a) shows AER model results, panel (b) 
AER/GSFC results, panel (c) AER/LLNL results, panel (d) GSFC results, and panel (e) LLNL results, 
all with PSC chemistry. 
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Figure 21 . Calculated perturbation in O 3 column (%) due to emission from HSCT aircraft in 2015 
relative to an atmosphere with only subsonic aircraft. Panel (a) shows AER model results, panel (b) 
AER/GSFC results, panel (c) AER/LLNL results, panel (d) GSFC results, and panel (e) LLNL results, 
all with PSC chemistry. 








